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c3 ■ Abstract 

We present a thermodynamic theory of plane coherent solid-sohd interfaces in multicomponent systems 

o 

en I subject to non-hydrostatic mechanical stresses. The interstitial and substitutional chemical components are 



treated separately using chemical potentials and diffusion potentials, respectively. All interface excess quan- 
tities are derived using Cahn's (1979) generalized excess method without resorting to geometric dividing 
surfaces. We present expressions for the interface free energy as an excess quantity and derive a general- 
ized adsorption equation and an interface Gibbs-Helmholtz equation which does not contain the interface 

^ ■ entropy. The interface stress tensor emerges naturally from the generalized adsorption equation as an ap- 

'^ . propriate excess over bulk stresses and is shown to be generally non-unique. Another interface property 

O I emerging from the generalized adsorption equation is the interface excess shear. This property is specific to 

coherent interfaces and represents the thermodynamic variable conjugate to the shear stress applied parallel 

'i> . to the interface. The theory reveals a number of Maxwell relations describing cross-effects between ther- 

'nI" I mal, chemical and mechanical responses of coherent interfaces. In Part II of this work this theory will be 

applied to atomistic computer simulations of grain boundaries. 

o 

m ■ PACS numbers: 64. lO.+h, 64.70.K-, 68.35.-p, 68.35.Md 
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I. INTRODUCTION 

Thermodynamics properties of interfaces can have a strong impact on microstructure de- 
velopment in materials by controlling phase nucleation, growth, coarsening and many other 
processes.-^. The interface thermodynamics developed by Gibbs^ was formulated in terms of 
interface excesses relative to an imaginary geometric dividing surface separating the coexisting 
phases. Gibbs defined the interface free energy 7 as the reversible work expended to create a 
unit area of the interface. He showed that, while other excess quantities generally depend on the 
choice of the dividing surface, 7 is unique and thus a meaningful and measurable quantity. Gibbs' 
work was focused on interfaces in fluid systems, for which he derived the adsorption equation 
expressing the differential ^7 in terms of differentials of temperature and chemical potentials of 
the components present in the system. 

Gibbs^ also discussed solid-fluid interfaces and pointed out that their interface area can change 
in two different ways: when a new area of the interface is formed at fixed thermodynamic states of 
the phases, and when the solid phase is elastically stretched parallel to the interface. The second 
process leads to the definition of the interface stress Tjj, a tensor quantity whose components are 
generally different from 7 and can be positive or negative.- Using a thought experiment with a 
solid equilibrated with three different fluids, Gibbs demonstrated that chemical potential of a non- 
hydrostatically stressed solid is not well-defined.^ At the time, solid solutions were unknown and 
Gibbs considered only single-component solids. When introducing 7 for solid-fluid interfaces, he 
specifically placed the dividing surface so that the interface excess of the solid component would 
vanish and there would be no need to talk about its chemical potential. As was recently pointed 
out,- this approach would not work for a multicomponent solid. 

Cahn^ proposed a more general form of the adsorption equation for hydrostatic systems by 
solving a system of Gibbs-Duhem equations for the bulk phases and for a layer containing the in- 
terface. By eliminating the Gibbsian construct of dividing surface, Cahn's method offers a greater 
freedom of choice of intensive variables in the adsorption equation. In particular, Cahn's formal- 
ism rigorously introduces the interface excess volume, a quantity which is by definition zero in 
Gibbs' thermodynamics. Cahn^ also proposed a Lagrangian (L) formulation of the Shuttleworth 
equation^ for phase boundaries, tlij = d'jL/deLij, and pointed to the importance of computing the 
derivative with respect to the elastic strain suj along a phase coexistence path. For solid-fluid inter- 
faces subject to non-hydrostatic mechanical stresses, Tjj has been formulated as an interface excess 



quantity^"— and computed by atomistic methods for several crystallographic orientations. ■^'^'i^iii 

In comparison with solid-fluid interfaces, thermodynamics of solid-solid interfaces is more 
challenging for at least two reasons. Firstly, such interfaces are capable of supporting shear stresses 
parallel to the interface plane. The interface response to such stresses depends on the degree of 
coherency and can vary from perfect sliding for fully incoherent interfaces to perfectly elastic re- 
sponse for fully coherent interfaces.— — The elastic response should obviously lead to additional 
terms in the adsorption equation, with coefficients representing what can be called "interface ex- 
cess shears". Such terms do not appear in existing formulations of interface thermodynamics.-^ 
Secondly, because of the undefined chemical potentials in non-hydrostatically stressed solids, a 
different treatment is needed for the — Fjfi/ij terms appearing in the adsorption equation for fluid 
systems.^ In the analysis of equilibrium between bulk solid phases, this problem was circumvented 
by using chemical potentials for interstitial components but diffusion potentials for substitutional 
components.— — It will be shown below that the same approach can be transferred to interface 
thermodynamic s . 

Elastic response of coherent interfaces was also analyzed within mechanical theories of 
interfaces.—^— In such theories, the interface is treated as a surface separating two elastic me- 
dia subject to applied stresses. Mechanical equilibrium conditions have been derived and possible 
excess deformations and stresses at both plane and curved interfaces have been identified. By 
contrast to thermodynamic theories,— — the mechanical theories do not impose the chemical or 
phase equilibrium conditions between the adjoining media. As a result, the interface free energy 
7 cannot be defined as the work of interface formation, which blocks the route to the formulation 
of the adsorption equation. 

In this paper we present a thermodynamic treatment of plane coherent solid-solid interfaces 
subject to mechanical stresses. Our analysis employs Cahn's generalized excess method^ and the 
solid-solid equilibrium theory developed by Robin^^ and Larche and Cahn.— — As the authors be- 
fore us,-^^— we do not rely on a priori chosen thermodynamic potentials; instead, all equations 
are derived directly from the First and Second Laws of thermodynamics. It is only after the deriva- 
tion is complete that some of the equations can be rewritten in simpler and/or more intuitive forms 
by introducing appropriate thermodynamic potentials. We start the paper by reviewing thermo- 
dynamics of a single solid phase (Section ^ and then formulate the coherent phase equilibrium 
conditions (Section HIH) in a form which prepares us for the subsequent thermodynamic analysis 
of interfaces. Section |IV] is central to this paper. In it, we define the free energy 7 of a coherent 



interface, reformulate it as an excess of appropriate thermodynamics potentials, introduce a num- 
ber of other interface excess quantities, and finally derive the generalized adsorption equation and 
the interface version of the Gibbs-Helmholtz equation. These equations identify and define the 
interface excess shear, a property specific to coherent interfaces and conjugate to the shear stress 
applied parallel to the interface. They also define the interface stress tensor as an excess quan- 
tity and demonstrate that it is not unique. We derive a number of Maxwell relations describing 
interesting cross-effects between different interface properties. In Section |V] we discuss how the 
proposed coherent interface theory can be applied to incoherent interfaces and grain boundaries. 
Finally, in Section |Vl] we summarize our work and outline possible future developments. 

n. THERMODYNAMICS OF A SOLID PHASE 

A. The network solid 

Our treatment of a solid phase is based on the concept of a network solid introduced by Robin^ 
and Larche and Cahn.— i^ We assume that the solid contains a penetrating network which is pre- 
served inside the solid and can be created or destroyed only at its boundaries. The network serves 
three functions: 

1. Enables a description of finite deformations of the solid'^ by associating physical points 
with network sites (or their small groups). 

2. Is capable of carrying mechanical loads, allowing the solid to reach mechanical equilibrium 
under non-hydrostatic conditions. 

3. Provides a conserved set of sites, called substitutional, which are completely, or almost 
completely, occupied by atoms. ^^ Accordingly, all chemical components can be divided 
into substitutional (residing on substitutional sites) and interstitial (otherwise).— 

In a crystalline solid, its lattice is formed by substitutional sites and satisfies all three network 
properties. Since our theory is intended primarily for applications to crystalline solids, we will 
adopt the terminology in which we refer to the network as "lattice" and to the network sites as 
"lattice sites". It should be noted, however, that our results are of more general validity and do not 
require that the solid have a long-range atomic ordering. The theory should be equally applicable 
to non-periodic structures such as network glasses or network polymers. 



B. Kinematics of deformation of a solid phase 

We will analyze the general case of finite deformations of a solid using the concept of a refer- 
ence state. ^^ The choice of the reference state is arbitrary, although it is often convenient to choose 
a stress-free state. We will use the same Cartesian coordinate system for both the reference and de- 
formed states. For any physical point defined by coordinates x[ (i = 1, 2, 3) in the reference state, 
its coordinates Xj in the deformed state are functions of the reference coordinates, x = x (x'). 
Any infinitesimal vector dxi connecting two physical points in the deformed state is related to 
the infinitesimal vector dx'j connecting the same two physical points in the reference state by the 
linear transformation 

dxi = ^ Fijdx'j, (1) 

i=l,2,3 

where tensor F is the deformation gradient with components 

^'^ dx'j ' 
It is assumed that J := det F 7^ and thus the reference coordinates can be expressed as functions 
of the deformed ones, x' = x' (x). The respective infinitesimal vectors are related by the inverse 
deformation gradient F^\ 

dx', = J2 P^i'dx,. (3) 

i=l,2,3 

Only six components of F are needed to completely describe all deformations (strains) of a 
solid. Without loss of generality, we will set all sub-diagonal components of F to zero. 



Fii F12 -F13 

F22 F2Z 

y F33y 



(4) 



with the Jacobian 

J = -F11F22-F33. (5) 

It can be shown that F~^ also has an upper-triangular form with diagonal elements 

Fr^ = l/Fu, ^ = 1,2,3. (6) 



The upper-triangular form of F implies that for any small volume element with the shape of a 
parallelepiped, its bottom and top faces remain normal to the X3 axis. Furthermore, the edge of the 
parallelepiped which is initially parallel to the Xi axis remains parallel to it during the deformation. 
Deformation of small volume element described by Eq. dH) is illustrated by a two-dimensional 
schematic in Fig.[TJ 

C. Thermodynamic description of a homogeneous solid phase 

Consider a homogeneous multicomponent solid containing K substitutional and L intersti- 
tial chemical components in a state of thermodynamic equilibrium. The defining property of the 
substitutional components is that their atoms fill a conserved set of lattice sites. Vacancies, i.e. 
unoccupied substitutional sites, are neglected for the time being and will be discussed separately 
in Section |Vll Thus, the total number N of substitutional atoms in any given reference region 
remains constant in all thermodynamic processes. Interstitial atoms occupy otherwise empty posi- 
tions between the lattice sites and their number in any given reference region can vary. Diffusion 
of both substitutional and interstitial atoms is allowed as long as it preserves the substitutional 
sites. 

Consider a homogeneous region of the solid containing a total of N substitutional and n inter- 
stitial atoms and obtained by elastic Odeformation of a homogeneous reference region of volume 
v. Suppose the reference region, and thus V and N, are fixed. Then the internal energy U of the 
region is a function of its entropy S, the amounts of individual chemical components Nk and ni, 
and the deformation gradient F: 

U = U{S,N,,...,NK,nu...,nL,F) (Fixed iV). (7) 

Due to the imposed substitutional constraint J2k ^k = N = const, only {K — 1) independent 
variations of N^ are possible. To implement this constraint, we can arbitrarily choose one of the 
substitutional components as the reference component and assume that each time we add to the 
solid an atom of a different substitutional component k, we simultaneously remove an atom of 
the reference component. ^^"^^ Let us choose component 1 as the reference and treat the amounts 
of all other substitutional components as independent variables. The amounts of the interstitial 
components can be varied without constraints. 

Consider a reversible variation of state of the solid, with a fixed reference region, when it 



exchanges heat with its environment, changes its chemical composition, and performs mechanical 
work. The differential of energy of the region is given by-^^^ 

K L 

fc=2 1=1 ij=l,2,3 

where T is temperature, fii are chemical potentials of the interstitial atoms and M^i are (K — 1) 
diffusion potentials of the substitutional atoms. According to Eq. ([8]), the diffusion potential M^i 
is the energy change when an atom of the substitutional component k is replaced by an atom of 
the reference component 1 while keeping all other variables fixed: 

dU dU 

In the last term in Eq. ([8]), P is the first Piola-Kirchhoff stress tensor, which is generally not 
symmetrical and is related to the symmetrical Cauchy stress tensor cr by^ 

P = JF-' ■ a- (10) 

(the dot denotes the inner product of tensors and vectors). Because F^^ is a upper-triangular 
matrix, the components P31, P32 and P33 are proportional to the respective components of cr: 

P-M = FnF22(yM = {J/F33) fT3^, ^ = l,2, 3, (11) 

where we used Eqs. (H), (|5]) and (|6]). 

In preparation for the analysis of interfaces in Section |W] we will rewrite the mechanical work 
term in Eq. ^ by separating the differentials dFu, dFi2 and dF22 from dFi^, dF23 and dF^^: 

K L 

dU = TdS + J2MkidNk + J2f^idni+ J] V FuF22(rs^dF,, + J] V'P,,dF,,. (12) 

fc=2 1=1 i=l,2,3 i,j=l,2 

The (K + L + 6) differentials in the right-hand side of this equation are independent and their 
number gives the total number of degrees of freedom of a homogeneous solid phase. 



D. Relevant thermodynamic potentials 

Various thermodynamic potentials can be derived from Eq. (fT2l) by Legendre transformations. 
As will become clear later, the potential relevant to coherent interfaces is 

K L 

^i:=U -TS -Y^MkiNk-Y^mi- Yl iVF,^/F^3)cT^i, (13) 

k=2 1=1 1=1,2,3 

where subscript 1 indicates the reference substitutional component and V = JV is the physical 
(deformed) volume of the homogeneous solid. For a fixed reference volume (and thus N), 

$1 = $i(r, M21, ..., Mki, /il, ..., /XL, f^Sl, ^^32, C^33, ^^11, ^12, ^^22)- (14) 

Using Eq. ([72)) we obtain 

K L 

ci$i = -SdT - Y, NkdMki - Y ^idf^i - Y (^^^3/^^33) da3^ + ^ V'Q^jdFj,, (15) 

fc=2 1=1 «=1,2,3 i,i=l,2 

where we denote 



Q := JF-i ■[^- Y y^^3n^I (16) 

\ m=l,2,3 ^^ / 

(I = 6ij is the identity tensor). Although Q is a 3 x 3 tensor, only its components Qu, Q21 and 
Q22 appear in Eq. (fTST ). 

While the potential $1 will prove to be useful in interface thermodynamics, its role in 
thermodynamics of a bulk phase is less obvious. Unless the state of stress is hydrostatic, 
this potential depends on the choice of the coordinate axes through the stress-strain variables 
^^315 cr325 cr335 -^11, -^12, -^22- In addition, $i depends on the choice of the reference state of strain. 

Eq. (fT3l) defines $i for a homogeneous solid region containing a given number N of substitu- 
tional sites. We can also define an intensive potential 0i as $i per substitutional atom: 

K L 

01 := $i/iV = U/N - TS/N - Y^hiCk - Yf'ici - ^ (^^*3/i^33) ^3^- (17) 

fc=2 1=1 i=l,2,3 

Here Ck := N^/N and q := rii/N are concentrations of substitutional and interstitial components 
per substitutional site, U/N, S/N and Vt are the energy, entropy and volume per substitutional site, 
respectively. 

Similarly, we can introduce K different potentials $„, and accordingly 0^, by choosing other 
substitutional components in as the reference species: 
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K L 



■.= <!>^/N = U/N-TS/N-J2MkmCk-J2fiiCi- J^ (^^,3/^33)^3^ (18) 



fc=l 1=1 1=1,2,3 



Note that we have extended the summation with respect to k from 1 to iC using the property M^k = 



0. Combining Eq. (1181) with known properties of diffusion potentials,— — namely Mik = —Mki 
and Mij = Mj^ + M^j, the following relationship between different 0-potentials can be derived: 

4>m - (t)n = Mmn, m,n=l,...K. (19) 

It also follows that 

K K K 

y^ MkmCk = ^ (0fc - 0m) C'fc = X] '^'^'^^ ~ '^™- ^^^^ 

fc=l fc=l fc=l 

Using Eqs. (fT8l ) and (|20l) . we obtain the following thermodynamic relation for a homogeneous 
non-hydrostatic solid phase: 

K L 

U-TS- Y, {VF^3/F33) ^3. = Yl "^^^^ + Yl /^'^'- (^1> 

j=l,2,3 fc=l «=1 

This equation closely resembles Gibbs' equation U — TS + pV = J^nf^nNn for hydrostatic 
systems (p being external pressure),- with (pk playing the role of chemical potentials. When the 
solid is in a hydrostatic state of stress, aij = —5ijp, the left-hand side of Eq. (|2TI) reduces to 
the Gibbs free energy U — TS + pV. Accordingly, 0^ become real chemical potentials of the 
substitutional components. 

E. The Gibbs-Duhem equation 

We can now derive a Gibbs-Duhem equation for a multicomponent non-hydrostatically stressed 
solid. To this end, we again consider a variation of state in which the solid region exchanges 
heat with its environment, performs mechanical work, and changes its chemical composition by 
switching chemical sorts of substitutional atoms (at fixed A^) and changing the amounts of inter- 
stitial atoms. Differentiating Eq. (fT3l) and using the relation d^i = Ndcpi and dU from Eq. (fT2l) . 
we obtain the following Gibbs-Duhem equation: 



K L 

0= -SdT-J2NkdMki-Nd<Pi-Y,nidfii- J^ {VF,^/ F^^) da^,+ ^ \/Uit/F,i. (22) 

k=2 1=1 1=1,2,3 «,j=l,2 

Applying Eq. (fT9l ). this equation can be rewritten as 

K L 

= -SdT - Y,Nkd(Pk - Y.'^idiii - J2 iVF,s/Fs3) da^, + J] V'Q,,dFj,. (23) 

A:=l i=l 1=1,2,3 i,i=l,2 

In the particular case of hydrostatic processes, Qij = while J2i=i 2 3 (^-^i3/-^33) d-^3i = 
—Vdp. In this case Eq. (|23l) reduces to the classical Gibbs-Duhem equation derived for fluids,- 

K+L 

= -SdT - J2 ^kdnk + Vdp, (24) 

fc=i 
where fXk are chemical potentials of chemical components and A^^ are their amounts. 

Eq. (|23l) is a relation between differentials of the intensive variables which characterize thermo- 
dynamic states of solids in equilibrium. By contrast to the standard Gibbs-Duhem equation (l24l) . 
it contains non-hydrostatic variations. 

ni. COHERENT EQUILIBRIUM BETWEEN SOLID PHASES 

A. Definition of coherency and coherent interface 

We next discuss coherent equilibrium between two homogeneous solid phases whose thermo- 
dynamic properties were introduced in Sectional] We assume that the two phases, which we refer 
to as a and (3, contain the same K substitutional and L interstitial components and are separated 
by an infinitely large, plane coherent interface normal to the X3 direction (Fig. |2]). Our definition 
of phase coherency follows the works of Robin^ and Larche and Cahn.— — Namely, a coherent 
transformation of a region of phase a to a region of phase (3 is accomplished by deformation of 
the lattice without creation or destruction of lattice sites. Thus, a coherent transformation fully 
preserves the reference region of the phase. All chemical components are allowed to diffuse dur- 
ing the transformation as long as the lattice sites remain intact. A more detailed discussion of the 
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concept of coherency and examples of coherent transformations can be found in Refs. 

If a coherent transformation occurs on one side of a plane selected inside a single-phase region, 
it produces a coherent interface between the old and new phases. Advancement of the transforma- 
tion front occurs by interface migration. For coherent phases, there is a single network of lattice 

10 



sites penetrating through both phases and deformed during the interface motion. In other words, 
the two-phase system can be described as a deformation map of the same reference region as 
each of the phases. Due to the lattice continuity across the interface, sliding is prohibited and the 
two-phase system responds to applied shear stresses elastically. (This is in contrast to incoherent 
interfaces, which do not support static shear stresses and do not preserve the lattice sites during 
their motion.) 

We will adopt the following kinematic description of coherent two-phase systems. The defor- 
mation gradients of the phases, F" and F^, are taken relative to the same reference state and have 
the upper- triangular forms. 



(25) 



F^ 
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-^13 





-^22 


-^23 


V " 
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F22 


F^ 

^23 


v» 





p/3 

-^33 y 



(26) 



where the superscripts indicate the phases. These forms ensure that the X3 direction in both phases 
remains normal to the interface plane during all deformations. In addition, the lateral deformation 
components Fn, F12 and F22 are common to both phases, which is a necessary condition for the 
absence of sliding. Thus, the two deformation gradients differ only in the components Fjs. The 
differences between these components form a vector. 



t : = 



p/3 _ -pa p/3 _ -pa pP _ pc 
-^13 -^13' -^23 -^23' -^33 -^33/' 



(27) 



which we call the transformation vector. Its geometric meaning is illustrated by the two- 
dimensional schematic in Fig.|2tc). 

B. Coherent phase equilibrium conditions 



The conditions of coherent phase equilibrium were derived for a single-component system by 
Robini^ and generalized to multicomponent systems containing both substitutional and intersti- 
tial atoms by Larche and Cahn^^^ (see Voorhees and Johnson^S for review). The equilibrium 
conditions can be summarized as follows: 
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(i) Temperature is uniform throughout the system. 

(ii) Diffusion potentials Mki of all substitutional components and chemical potentials fii of all 
interstitial components are uniform throughout the system. 

(iii) The internal mechanical equilibrium condition, V ■ P = 0, is satisfied inside each phase 
(the divergence is taken with respect to the reference coordinates). 

(iv) The traction vector is continuous across the interface, 

n'" ■ P'^ = -n'" ■ P^ (28) 

where vectors n'"" and n = — n'" are unit normals to the phases in the reference state. ^^ This 
condition reflects the continuity of the displacement vector across the interface. From Eqs. (|28] 
and (fTT)) it follows that the Cauchy stress components aai, (732 and ass are also continuous across 
the interface. 

(v) Finally, the so-called phase-change equilibrium condition^^^ must be satisfied. This con- 
dition expresses equilibrium with respect to virtual displacements of the interface in which a layer 
of one phase reversibly transforms to a layer of the other. Rewritten in our notations, the phase 
change equilibrium condition derived by Larche and Cahn^^^ reads 

K L 



k=2 1=1 i=l,2,3 

K L 



k=2 1=1 1=1,2,3 

Here U, S, V are the energy, entropy and volume of the phases obtained by deformation of the 
same reference region. The total number of substitutional atoms is equal in both phases, A^" = A^'^, 
whereas the total number of interstitial components can be different (n" ^ n^). 

The equilibrium conditions (i)-(iii) are common to all types of interfaces. The differences 
between the coherent, incoherent and other types of interfaces lie in the remaining conditions (iv) 
and (v). 

C. Derivation of the phase-change equilibrium condition 

The phase-change equilibrium condition (|29l ) was obtained from equation (41) of Larche and 
Cahn^"^ by inserting our upper-triangular deformation gradients ([25]) and (|26l) and the interface 
normal n = (0, 0, 1). Note that Eq. ^ contains the terms (^^-Ffs/^33 - ^"^13/^33) ^31 and 
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^''-^^/-^ss ~ ^"-^23/-^33 ) "^32 proportional to the shear stresses 0-31 and (T32. These terms are 
specific to coherent interfaces and vanish for incoherent, solid-fluid and fluid-fluid systems which 
do not support such stresses. To elucidate the meaning of these terms and set the stage for the 
analysis of interface thermodynamics, we will present an alternate derivation of Eq. (1291 ) which 
assumes that the equilibrium conditions (i) through (iv) are already satisfied. 

At fixed values of the intensive variables T, M21, ..., Mki, fii, ..., /iL, 0-31, 0-32, (733, Fn, F12, -F22,— 
equilibrium between the phases is neutral, i.e., the interface can reversibly migrate up and down 
without altering thermodynamic states of the bulk phases. The phase change equilibrium 
condition expresses the neutrality of this equilibrium with respect to interface displacements. 
Consider a homogeneous layer of phase a parallel to the interface and containing a total of 
A^ substitutional atoms. Suppose the interface traveling down passes through this layer and 
transforms it completely to a layer of phase /3. The initial and transformed states of the layer are 
shown schematically in Fig. I2a,b). In both states, the layer contains the same total number of 
substitutional atoms, whereas the total number of interstitial atoms can be different. 

Let us compute the change in internal energy of this layer. Because the transformation is re- 
versible, this change depends only on the initial and final states (i.e., homogeneous phases a and 
/3) and not on the transformation path. As the interface traverses the layer, it creates intermediate 
states that are not homogeneous. Instead of examining this actual transformation process, we will 
consider another, imaginary path on which the transformation occurs by homogeneous deforma- 
tion of the layer with a simultaneous change in its chemical composition. Since the layer remains 
homogeneous during this process, its energy change can be obtained by integrating Eq. (fT2] ) de- 
rived previously for homogeneous variations. Remembering that the intensive parameters are 
fixed, the integration gives 

K L 

k=2 1=1 (30) 

+ J2 (y^Ffs/Fis - y^FtjFi^ ^3.. 

j=l,2,3 

The last term in Eq. (fT2l) does not contribute to this equation because Fn, F12 and F22 are not 
varied. Eq. (|30l) recovers the phase-change equilibrium condition (|29l ). 

This derivation emphasizes that the last term in Eq. (1301 ) represents the mechanical work, Wm, 
done by the stress components 0-34 during the phase transformation. This work term can be rewrit- 
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ten as 



W, 



m 

1=1,2,3 



J] F11F22V' (i^5 - Fi) asi = FuF22V'a- ■ t, (31) 



where t is the transformation vector defined by Eq. (|27t and illustrated in Fig.[2tc). It is important 
to note that, while V^F^/F^^ and V^FJ^/F^r^ individually depend on the choice of the reference 
state of strain, vector t is an invariant and in principle measurable quantity characterizing the ge- 
ometry of the transformation.^'^ For incoherent and other interfaces incapable of supporting static 
shear stresses, Wm reduces to -Fii-F22^'c"33i^3 = {V^ — V") o"33- For coherent interfaces, addi- 
tional work is done by the shear stresses along the components of t projected on the interface 
plane. 

Using the thermodynamic potential 0i defined by Eq. ([17]) . the phase-change equilibrium con- 
dition (|29l) can be formulated as simply (f)f = (f)^. Furthermore, by choosing other substitutional 
components as reference species, the following K relations can be obtained: 

C = 0^:=0m, m = l,...,K. (32) 

Thus, in a system with K substitutional chemical components, there are K potentials that have the 
same value in coexisting phases. This result resembles Gibbs' condition of equilibrium between 
fluid phases,^ with (pm playing the role of chemical potentials. 

D. The equation of coherent phase coexistence in the parameter space 

The Gibbs-Duhem equation (|22)) establishes a relation between the differentials of (K + L + 7) 
intensive parameters characterizing a single-phase solid under stress. When two solid phases 
coexist, their equilibrium imposes an additional constraint on possible variations of state of the 
phases. This constraint can be formulated by writing down the Gibbs-Duhem equation for each 
phase in terms of the same set of intensive parameters and requiring that the two equations hold 
simultaneously: 

K L 

fc=2 1=1 (33) 

- E (^"^3/^3^3) da,. + E V'-Qf.dF,,, 

i=l,2,3 «,i=l,2 
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K L 

= -Sf^dT - J^^kdMki - N^d(t)i - ^nfd^i 

k=2 1=1 (34) 

j=l,2,3 i,i=l,2 

Note that these equations are written for arbitrarily chosen amounts of the phases, i.e. generally 
iV" ^ N^. They can be combined into one equation by eliminating one of the differentials. This 
elimination leads to the equation 

K L 

= -{S}xdT - Y,{Nk}xdM,, - {N}xd<P, - Y.Mxdf^i 

k=2 1=1 (35) 

- J2 {VF,s/F,s}xdas^ + J2 {V'Q^,}xdF^^, 

1=1,2,3 «J=1,2 

where X is one of the extensive properties S, N^ (k = 2, ..., K), N, ni (I = 1, ..., L), VFi^/F^^ 
(i = 1, 2, 3) or V'Qij (i,j = 1,2). The curly braces are defined by 

{Z}x ■■= ^" - Z^X"/X^ (36) 

for any pair of extensive properties Z and X. The physical meaning of {Z}x is the difference 
between the property Z of the two phases when they contain the same amount of X. For example, 
{S}n is the difference between entropies of two homogeneous regions of the phases containing 
the same total number of substitutional atoms. 

For any choice of X out of the above list, the respective differential coefficient in Eq. (|35l) van- 
ishes because {X}x = 0. The remaining (K + L + 6) terms form a differential equation defining 
the coherent phase coexistence hypersurface in the configuration space of intensive parameters. 
Thus, a system of two coexisting coherent phases is capable of (i^ + L + 5) independent vari- 
ations, which is one degree of freedom less than for each phase taken separately. Knowing one 
equilibrium state of the two-phase system, all other states can be found by integrating Eq. (|35]) 
along different paths on the phase coexistence hypersurface. 

Eq. (|35l) is an important result of this paper. It provides the phase rule for equilibrium between 
coherent phases and offers flexibility in choosing the independent variables corresponding to the 
available degrees of freedom through the choice of X. It generalizes the equation of phase coexis- 
tence derived by Gibbs for solid-fluid interfaces^ by incorporating shears parallel to the interface. 
Such shears are represented by the additional terms {V^-Fjs/Fssjxc^crsj with i = 1,2. To further 
elucidate the physical meaning of these terms, consider coherent equilibrium between two binary 
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substitutional solid solutions. For variations of the shear stress at a constant temperature and fixed 
lateral dimensions of the system, 

dM2i ^ {VFjs/FssJN , 

dasi {NkJN 

where we chose X = N. This relation predicts that to maintain the equilibrium, variations in 
the diffusion potential in response to variations in the shear stress must be proportional to the 
transformation shear and inversely proportional to the difference between the phase compositions. 
In other words, this relation describes changes in the phase compositions caused by applied shear 
stresses. 

It should be emphasized that Eq. (|35]) has been derived under the assumption of interface co- 
herency. One might think that the phase coexistence equation for incoherent interfaces could be 
obtained as simply a particular case of Eq. (|35T ) when a^i and 0-32 are zero. This is not so. In the ab- 



sence of coherency, the lateral deformations of the phases (F^ and F[-,i,j = 1 , 2) are not required 
to be equal and can be varied independently. For example, one of the phases can be stretched in 
a certain direction parallel to the incoherent interface while the other compressed in the opposite 
direction. This deformation produces interface sliding, which is a possible process for incoherent 
interfaces. Furthermore, because the lattice sites can now be created or destroyed when one phase 
transforms to the other, the deformation gradients F° and F^ must be defined relative to dijf event 
reference states. The incoherent phase coexistence equation would have to be re-derived from the 
start, which is beyond the scope of this paper. 

IV, INTERFACE THERMODYNAMICS 

A. The interface free energy 7 

We are now ready to analyze thermodynamics of coherent interfaces. In this section we derive 
expressions for the interface free energy 7 defined as the reversible work expended for creation of 
a unit interface area. As above, we imagine two coexisting phases a and (3 separated by a coherent 
plane interface (Fig. |2l), but we now include the interface region as part of the system. Recall that 
the deformation gradients F" and F^ were previously introduced for homogeneous phases and 
remain undefined within the highly inhomogeneous interface region. We therefore need to devise 
a method for introducing 7 and other interface excess quantities without defining a deformation 
gradient inside the interface region. 
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As discussed earlier, the coherent two-phase equilibrium is neutral when the intensive param- 
eters T, M21, ..., Mki, Hi, ..., fJ'L, 0-31, 0-32, 0-33, Fii, F12, F22 are fixed.— Consider a homogeneous 
region of phase a in the shape of a parallelepiped with a reference volume V. Two faces of the 
parallelepiped are parallel to the interface and one edge is parallel to the xi axis. This is illus- 
trated by a two-dimensional schematic in Fig. Sb), where the parallelepiped is represented by a 
parallelogram. Suppose the interface spontaneously migrates and enters this region, turning it into 
an equilibrium two-phase system [Fig. HJc)]. Due to the coherency condition, the cross-section 
of the region parallel to the interface remains the same at every height. However, the shape of 
the region changes due to the phase transformation strain. Consider a particular position of the 
interface inside the two-phase region such that the upper and lower boundaries of the region are 
deep inside the homogeneous phases not perturbed by the presence of the interface. Suppose the 
lower boundary of the region is fixed. Then the position of the upper boundary generally changes 
as as result of the phase transformation. Denote the displacement vector of the upper boundary B. 

The geometric meaning of vector B is illustrated by the two-dimensional schematic in Fig. |3] 
The initial region abed is a deformation map of a reference region a'h'c'd! with the deformation 
gradient F" [Fig. |3]^a,b)]. After the upper part of the region transforms to phase (3, it becomes 
a map of the corresponding upper part of the reference region with the deformation gradient F'^. 
The reference comers c' and d! are thus mapped to some physical points c* and d* within the (3 
phase [Fig.|3tc)]. Vector B is defined as cc*, or equivalently, dd* [Fig.jS^d)]. Note that due to the 
conservation of sites by coherent interfaces, the two-phase region contains the same number A^ of 
substitutional sites as the initial region of phase a. 

Vector B is used for calculation of the mechanical work Wm performed by stresses when the 
discussed region transforms to the two-phase state. Since the cross-section of the region remains 
fixed, the mechanical work is done only by the stress components a^i when the upper boundary 
is displaced by vector B. Thus, Wm = An°'-(T ■ B, where n° is the unit normal to the interface 
pointing into phase /3, n"-cr is the traction vector, and A is the cross-sectional area. 

To keep similarity with the mechanical work terms derived previously for homogeneous phases 
(e.g., Eq. (|30l) ). we want to express Wm through some deformation gradient. To this end, we 
formally define a homogeneous deformation gradient F relative to the same reference state as used 



17 



for the homogeneous phases: 

^ Fn Fi2 {F,% + B,A'/V')^ 
F:= F22 {F^^ + B2A'/V') , (38) 

^ iF,% + B,A'/V')^ 

where A' is the cross-sectional area of the interface in the reference state. We will refer to F 
as the "average" deformation gradient of the region. The geometric meaning of F is the affine 
transformation that carries the parallelepiped representing the reference region of phase a to the 
parallelepiped formed by the comers of the two-phase region after the phase transformation. In 
the two-dimensional schematic shown in Fig. [3l F transforms the reference region a'b'c'd' to the 
parallelogram abc*d*. The latter is shown separately in Fig. IHe). It should be noted that both B 
and F generally depend on the choice of the reference thickness V /A' of the a phase region and 
on the position of the interface within the two-phase region. In terms of F, the mechanical work 
term can now be rewritten as 

W^ = AJ2 ^3^B^ = J2 (y^^s/F^s ' V^F^JF^^ asi, (39) 

1=1,2,3 1=1,2,3 

where V"' = FuF22F^^V' and V = -^11-^22-^33^^' are physical volumes of the a phase region and 
the two-phase region, respectively.— 

We next calculate the change in internal energy of the region when it reversibly transforms 
from phase a to the two-phase state [Fig. [3];b,c)]. Instead of tracking the actual motion of the 
interface into the region, we will consider only the initial and final states and imagine a differ- 
ent reversible process between them. Specifically, consider a process of homogeneous phase 
transformation a — )■ /3 in the upper part of the region at fixed N and fixed intensive parame- 
ters T, M21, ..., Mki, /ii, ..., /iL, C31, o"32, o"33, -^11, -^"12, -^22- Sincc the transformation occurs in an 
open system, its energy changes due to the following processes: (i) heat exchange with the envi- 
ronment, (ii) diffusion of atoms in and out of the system (at constant A^), (iii) mechanical work 
Wm performed by stresses applied to the boundaries of the region, and (iv) non-mechanical work 
Wnm associated with local atomic rearrangements leading to the formation of the interface. Using 
Eq. dill) for Wm, we have 
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K L 

k=2 1=1 (40) 

i=l,2,3 

where the extensive quantities with and without superscript a refer to the initial and final states, 
respectively. 

We define the interface free energy 7 as the non-mechanical work done per unit interface area, 
i.e., 'J A := Wnm- Using Eq. (fTT] ) for the a phase, Eq. (|40l ) can be simplified to 

K L 

jA = U-TS-J2MkiNk-<piN-J2f^ini- J] {VF.s/Fss) a3^, (41) 

fc=2 1=1 1=1,2,3 

or expressing the diffusion potentials through the 0-potentials using Eq. (fT9l) . 

K L 

'yA = U-TS~Y.'^kNk-Y,f'ini~ Yl (^^..3/^33)^3.. (42) 

k=l 1=1 i=l,2,3 

These equations can be rewritten in a shorter form by introducing the $1 potential of a two-phase 
region by analogy with Eq. ([T3] ): 

K L 

<l>i:=f/-T5- J^MfciiVfc- J]/i,n,- J] (^^^,3/^33)^3,. (43) 

k=2 1=1 i=l,2,3 

Then 

^A = ^i-N(l)i, (44) 

so that 7 is an excess of the $1 potential per unit interface area. Of course, instead of component 
1 we could have chosen any other substitutional component as a reference. 

Eqs. (|4TI) and (l42l) express the total interface free energy 7^4 through properties of an arbitrary 
region containing the interface. While 7A is uniquely defined by these equations, the individual 
terms appearing in the right-hand side depend on the location of the boundaries of the region. To 
express these terms through interface excesses that are independent of the boundaries, we need to 
subtract the contributions of the homogeneous phases. To this end, we select two arbitrary regions 
inside the homogeneous phases. Such single-phase regions can be chosen either inside or outside 
the two-phase region. The latter case is illustrated in Fig. |4l Let the total numbers of substitutional 
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atoms in the single-phase regions be A^" and A^'^, respectively (generally, N""^ N^). Eq. ([TT 
applied to these regions gives 



K 



= f/- - TS" - ^MfciiV" - 0iiV" - ^/i,< - Y, {VF,,/F,3ra,i 



k=2 



1=1 



i=l,2,3 



and 



(45) 



K L 

= U^ -TS^ -J^^f'^^k -<PiN^ -J2f^^i ' E {VF,s/F,sfas^. (46) 

k=2 1=1 i=l,2,3 

Eqs. PT]) . (|45l ) and (|46l ) form a system of three linear equations with respect to the same intensive 
variables. We solve this system of equations for 'jA using Cramer's rule of linear algebra.- The 
solution has the form 



K 



jA = [U]xY - T[S]xY - J2Mki[Nk]xY - MN]xY - J]/i/K]xy - Yl [^^^s/i^aslxy^s^, 

fe=2 1=1 1=1,2,3 

(47) 
where X and F 7^ X are any two of the extensive quantities U, S, N^ (k = 2, ..., K), N, ni 
(/ = 1, ..., L) or VFis/F^s {i = 1,2, 3). Note that the last member of this list, corresponding to 
i = 3, is simply volume V. Using Eq. (fT9l) we obtain the equivalent form of 'jA: 



K 



7A = [U]xY - T[S]xY - 5^0fc[iVdxy - 5^/Wzh]xy - Y. [^^*3/i^33]xya3.. (48) 



k=l 



1=1 



i=l,2,3 



The coefficients [Z]xy are computed as ratios of two determinants:- 



[Z] 



z 


X 


Y 


Z" 


x^ 


yet 


z^ 


x^ 


Y^ 



XY ■' 



(49) 



X^ Y^ 

The quantities in the first row of the numerator are computed for the region containing the inter- 
face, whereas all other quantities are computed for arbitrary homogeneous regions of phases a and 
(3. By properties of determinants. 
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[X]xY = [Y]xY = 0, (50) 

so that two terms in each of the Eqs. (|47l) and (l48l) automatically vanish. 

The coefficient \Z]xy has the meaning of the interface excess of extensive property Z when 
the region containing the interface contains the same amounts of X and Y as the two single-phase 
regions combined; in other words, when the excesses of X and Y are zero. Thus, the excess of any 
property Z is not unique; it generally depends on the choice of the reference properties X and Y . 
If either X or F is volume, then \Z]xy has the meaning of the excess of Z relative to a dividing 
surface similar to Gibbs' formulation of interface thermodynamics.- 

The excesses [A^fejxy and [n/]xy characterize the segregated amounts of substitutional and 
interstitial components, respectively. The terms [V"]xy, [^-^isZ-^sslxy and \VF2z/ Fyi\xY define 
the excess volume and two excess shears, respectively. For example, 

yF,3/F33 s V 

{yF,,/F^,f s^ v^ 



[FFis/Fsa] 



sv 



1,2. 



(51) 



In this case, the excess shears are taken with respect to a dividing surface for which the excess 
of entropy is zero. The excess shears are properties specific to coherent interfaces. They have 
no significance for incoherent solid- solid interfaces, solid-fluid interfaces or any other interfaces 
that cannot be equilibrated under applied shear stresses. By contrast, the excess volume [\^]xy is 
common to all types of interfaces.- The numerical values of the excess volume and excess shears 
depend on the choice of the reference properties X and Y. In the Gibbian formalism of dividing 
surface the excess volume is zero by definition. 

The total interface free energy 'jA can be expressed through excesses of different thermody- 
namic potentials corresponding to possible choices of X and Y. As already noted, ^A can be 
expressed as an excess of potential $i. Using our square bracket notation. 



K L 

jA=[U-TS- Y^^kiNk - 5^/^;nz - J^ (^^,3/^33) au]NV = ['^i]nv, (52) 

fc=2 1=1 i=l,2 

i.e., the excess of $1 must be taken relative to the dividing surface for which the excess of the total 
number of substitutional atoms is zero. As another example. 
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K L 

7/1 = [t/ - J]MfciiVfc - ^/i^n^ - Y, {V^^3/F33) (Ts,]ns, (53) 

k=2 1=1 i=l,2,3 

i.e., 'J A is an excesses of the potential appearing in the square brackets when the excesses of the 
total number of substitutional atoms and entropy are zero. The flexibility in expressing the same 
quantity 'jA through excesses of different thermodynamic potentials can be useful in applications 
of this formalism to experimental measurements and simulations. 

B. The adsorption equation 

Having introduced the interface free energy, we are now in a position to derive the generalized 
adsorption equation for coherent interfaces. As the first step, we will compute the energy differ- 
ential dU for a two-phase region containing the interface. We will take a region in the shape of a 
parallelepiped as shown schematically in Fig.jSj^e). Recall that this shape is a map of the reference 
region of phase a containing the same number of substitutional atoms as in the parallelepiped. 
This deformation map is formally defined by the deformation gradient F given by Eq. (|38l) . Con- 
sider a reversible variation in which this region exchanges heat and atoms with its environment (at 
fixed A^) and performs mechanical work by elastically changing its shape and dimensions. The 
mechanical work, dWm, is done by the stresses applied to all faces of the parallelepiped and equals 
the sum of the total forces exerted on the faces times their displacements. The calculations give 

dWm = Yl y'FnF22cr3^dF,, + J] V'PijdF^i, (54) 

i=l,2,3 j,i=l,2 

where P := JF ■ ct is a formal analog of the first Piola-Kirchhoff stress tensor, W is the true 
stress tensor averaged over the volume of the parallelepiped, J := det F, and V is the reference 
volume of the phase a region. Because the stress components a^i are coordinate-independent and 
the lateral stress components ctjj (i, j = 1, 2) depend only on the coordinate X3, it is only the lateral 
stress components that must be averaged over 0:3 in order to obtain W. Using the above expression 
for dWm, the energy differential equals 

K L 

dU = TdS + Y^kidNk + Yl^i^^i 

k=2 1=1 (55) 

+ Y y'FuF22cr3dFi3 + Y V'PrjdFji. 

i=l,2,3 «J=1,2 
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This equation looks similar to the previously derived Eq. ((121) and constitutes its generalization to 
inhomogeneous systems containing a coherent interface. 

At the next step, we take the differential of Eq. (|4TI) and insert dU from Eq. (|55] ). After some 
rearrangement we obtain 

K L 

d{jA) = -SdT -^NkdMki - Ndcpi -Y^mdfii 

k=2 1=1 (56) 

1=1,2,3 i,i=l,2 

where we introduced the tensor 

Q:=JF".(a- J2 f^^3„.l). (57) 

\ m=l,2,3 33 J 

For a homogeneous phase, Q reduces to the earlier introduced tensor Q, see Eq. (fT6l) . 

The differentials in the right-hand side of Eq. (1561 ) are not independent. There are two con- 
straints imposed by the Gibbs-Duhem equations (l33l) and (|34|) containing the same differentials. 
Solving the system of equations (l56l) . (|33l) and (|34l) by Cramer's rule, we finally obtain the gener- 
alized adsorption equation 

K L 

di'yA) = -[S]xYdT-J2[Nk]xYdMk,-[N]xYd<j),-Y.i''^^^y'^^i 

k=2 1=1 (58) 

- J2 [yF^3/F33]xYda3^ + J^ [V%]xYdF,i, 
1=1,2,3 i,j=l,2 

where X and F are two of the extensive properties 5*, Nf^ (k = 2, ..., K), N, rii {I = 1, ..., L), 
VF,,/F,s (^ = 1, 2, 3) or V% {t,j = 1,2). 

Note the significant difference between Eqs. (l56l) and (1581) written for the same differential 
d{'yA). In Eq. (|56l) . the differential coefficients are properties of the entire region containing 
the interface. These properties depend on the choice of the boundaries of the region and thus 
have no physical significance. In the adsorption equation (l58l) . on the other hand, the differential 
coefficients are interface excesses [Z]xy defined by Eq. (|49l) . For a given choice of the reference 
properties X and Y, such excesses are independent of the boundaries of the region. Furthermore, 
the number of differentials in the right-hand side of Eq. (l56l) exceeds the number (K + L + 5) 
of degrees of freedom of a coherent two-phase system predicted by Eq. (l35l) . By contrast, due 
to the property (|50] ) of determinants, two terms in Eq. (|58) automatically vanish, leaving exactly 
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(K + L + 5) independent differentials. Each of the remaining excesses [Z]xy can be expressed 
as a partial derivative of 'jA with respect to the corresponding intensive variable and is therefore a 
measurable physical quantity. In terms of the 0-potentials, the adsorption equation takes the form 

K L 

d {'J A) = -[S]xYdT - "^[Nklxydcpk - ^K]xyc?/i/ 

_ _fc=i 1=1 _ (59) 

- Yl [yF,s/F,3]xYdas, + J^ [y^UydF,,. 

4=1,2,3 ij = l,2 

The adsorption equation corresponding to Gibbs' formalism of the dividing surface is obtained 
as a particular case of our adsorption equation when either X = V orY = V. Although the excess 
volume [V]xY disappears, the excess shears [VFi^/F^^lxY and [I^-Fjs/Fssjxy still remain. These 
additional terms are not present in Gibbs' interface thermodynamics^ or in Cahn's work.- 

C. The interface stress 

The terms in the adsorption equation that contain differentials of the lateral deformation com- 
ponents Fii, Fi2 and F22 represent contributions to 7A coming from elastic deformations of the 
interface. These terms define the interface stress, the quantity which was first discussed by Gibbs 
in the context of solid-fluid interfaces .- 

To formally define the interface stress tensor, choose the current state of one of the phases as 
the reference state of strain. Then Fn = F22 = 1, -^12 = 0, and Eq. (l58l) becomes 

K L 

d (7A) = -[S]xYdT - Y,[Nk]xYdMki - [N]xYd(t)i - J^Nxyrf/x; 

k=2 1=1 (60) 

- Y. [yFr^/F^Yda:,, + ^ '^^^ ^de,,, 

4=1,2,3 i,i=l,2 

where 

r^/ ■■= \[V'Qn]xY, T,Y := ^[^'^22]^^, rj^^ = r,Y := j[V'Qu]xy. (61) 

In Eq. (l60l) . de is a (2 x 2) small-strain tensor with components den = dFu, de22 = dF22 and 
dei2 = de2i = ^dFi2. The (2 x 2) symmetrical tensor r defined by Eqs. (|6T] ) is the interface stress 
tensor describing changes in the interface free energy due to its elastic deformations. 

As other interface excess quantities, r generally depends on the choice of the reference prop- 
erties X and Y and is therefore not unique. However, in the particular case when both phases are 
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in a hydrostatic state of stress under a pressure p, r becomes independent of X and Y and is given 
by 



'^ij = -J (^ij + P^^j) ■ ^, J = 1, 2. (62) 



V 

'a 

This equation immediately follows from the definition of the square bracket [V'Q^j]xy and the 
fact that for hydrostatic phases Q^j = Q^ = 0. 

Eqs. (|6TI) provide a recipe for interface stress calculation when the phases are subject to non- 
hydrostatic stresses, particularly when such stresses are different in the two phases (e.g., when one 
phase is under lateral tension while the other under lateral compression). Previous calculations 
of interface stresses were focused on unstressed or hydrostatically stressed phases. For solid- 
fluid interfaces, the calculations for hydrostatic phases employed equations similar to (|62] ). Non- 
hydrostatic stresses were included only in surface stress calculations in single-phase systems.- 
Calculations of r between non-hydrostatic solid phases using Eqs. (1611 ) is an uncharted territory 
and could be addressed in future work. 

In the remainder of this paper, the lateral deformations of a two-phase system will be described 
by the small-strain tensor de instead of the lateral components of the deformation gradient. As 
mentioned above, this implies that the current state of one of the phases is chosen as the reference 
state of strain. It should be emphasized that (i) this assumption only reflects a particular choice 
of the kinematic description of deformations, not a physical approximation, and (ii) the normal 
and shear components Fj3 describing the transformation strain between the two phases can still be 
finite. 

D. Lagrangian and physical forms of the adsorption equation 

Until this point we dealt with total excess quantities related to the entire interface with an area 
A. It is often useful to define specific excesses, i.e. excesses per unit interface area in either the 
current state or the reference state of strain. In the former case, the excess quantity is referred to as 
physical while in the latter case as Lagrangian. For example, 7 is the physical specific excess of the 
interface free energy. The interface stress defined by Eq. (|6TI) is the physical specific excess of the 
tensor quantity V'Q^y One can also define the Lagrangian interface free energy, 7/, := (7 A) /A', 
and the Lagrangian interface stress, 

^Lii '■= -^,[y'Qii\XY, rf^ ■■= -Jj[V'Q22]xY, T^u = ^l-K '■= ^[^'Qi2]xy- (63) 
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rS=[^\ . (64) 



6 



Using Eq. (1631) and the adsorption equation (1581) . we obtain the relation 

XY 
L \ 

'^^^- " \de 
A similar relation involving the Lagrangian 7 was proposed by Cahn.^ Here, the superscript XY 
in the right-hand side indicates that the partial derivative is taken at fixed intensive parameters, 
other than Cij, that appear in the adsorption equation when X and Y are chosen as the reference 
properties. 

The Lagrangian form of the adsorption equation is obtained by dividing Eq. (|60l ) by A': 

K L 

, [o\XY ,^ sr^[^k\xY ,,. W\xY ,, sr^[ni]xY, 

k=2 1=1 (65) 



_k=2 

E\VFi'i/F^^]xY J . sr^ XY J 
da-i^+ 2^T^^deji 



A 

i=l,2,3 «J=1,2 



where the differential coefficients are Lagrangian specific excesses. The physical form of the 
adsorption equation is obtained by differentiating 'jA in Eq. (|60l ) and using the relation dA = 

^ z2i,j=i,2 ^ijdei 



-"i-r 



d^ = J-^dT - f^^-^dM,, - i^rf0, - x: ^^d^^ 



A ^ A ^^ A ^' ^^ A 

_fc=2 1=1 

E\VFi'i/F^^]xY , . yr^ / xy x ^^ 
J dcr3^+ 2^ [T.y - dijj) dCji. 

i=l,2,3 ij=l,2 



(66) 



Now the differential coefficients give physical specific excesses. From this equation, we immedi- 
ately obtain the generalized form of the Shuttleworth equation:- 

The original Shuttleworth equation"^ was derived for an open surface of a stress-free single- 
component solid deformed isothermally. Eq. (|67l ) has been derived for coherent interfaces in 
multicomponent systems in an arbitrary state of stress. It actually represents a set of equations 
corresponding to different choices of X and Y and thus different deformation paths. 

Just as the Shuttleworth equation describes the effect of lateral strains on the interface free 
energy, the following equations describe the effect of the shear and normal stresses of 7 and 72, in 
the physical and Lagrangian forms, respectively: 



XY 



djY' [VF,,/F,3]XY 



dasi J A 
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i = 1,2,3, (68) 



^d^J ^ 1^ ' ^ = ^'2' 2- ^^^^ 

E. Thermodynamic integration 

We will now derive another version of the adsorption equation that can be useful in applications. 
In principle, the interface free energy 7 can be computed by integration of the adsorption equation 
along a phase coexistence path knowing an initial value. However, the excess entropy \S\xy 
appearing in this equation is rarely accessible by experiments or simulations. To avoid calculation 
of [•S'Jxy, we can eliminate it by combining Eqs. (l47l) and (l60l) to obtain 



ib^\ . _ffl^,r-5:i%.M..-M^.^,-5:M?r,,, 



fc=2 l=\ 



;=1,2,3 i,j=l,2 



(70) 



where the thermodynamic potential ^ is defined by 

K L 

$ := f/ - ^iV.Mfci - iV0i - Y^^id^ii - Y, (^3^VF,s/F33. (71) 

fe=2 1=1 2=1,2,3 

It is straightforward to derive physical and Lagrangian forms of this equation, whose left-hand 
sides will contain d (7/T) and d ('Jl/T), respectively. 

In the particular case when only temperature is varied, Eq. (TTOl) gives 

XY 

dT J T2 

This equation is similar to the classical Gibbs-Helmholtz equation^ 



dijA/Tiy ^J3^_ (72) 



d{G/T) \ U + pV 

dT )^ T2 ^ ^ 

for single-component fluid systems. Eq. (fTOl) can be viewed as a generalization of the Gibbs- 
Helmholtz equation to interfaces in multicomponent systems. 

Eq. (ITOI) can be used to compute 7 by integration of 'yA/T along a trajectory on the phase 
coexistence hypersurface in the configuration space of variables. The advantage of this integration 
is that it does not require knowledge of [SJxy • Free energy calculations for solid-solid interfaces 
in multi-component elastically stressed systems are presently non-existent and could be initiated 
by applying the proposed thermodynamic integration approach. 
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F. Maxwell relations 

Because the adsorption equation contains the perfect differential of 7A, it generates a number 
of Maxwell relations between partial derivatives of the excess quantities. Similarly, the Gibbs- 
Helmholtz equation (TTOl ) is the perfect differential of 7A/T and also generates Maxwell relations. 
We will focus on Maxwell relations that involve the effects of mechanical stresses and strains 
on interface properties. For hydrostatic precesses, such relations were discussed by Cahn.^ The 
additional terms in the adsorption equation introduced in this work, such as the variations in the 
shear stresses (731 and (732, lead to a number of additional Maxwell relations. The Lagrangian 
and physical forms of the adsorption equation produce different Maxwell relations, which will be 
presented below side by side. In the partial derivatives appearing in these relations, the variables 
which are held constant are dictated by the particular choice of the extensive variables X and 
Y. Thus, each Maxwell relation actually represents a set of relations corresponding to different 
choices of X and Y. 

1. Mechanical relations 

The first set of Maxwell relations examines how the lateral deformations dcij and the stresses 
a^k affect the interface excess volume, excess shears and interface stress. Using the Lagrangian 
and physical forms of the adsorption equation, we obtain 



dr^ _ dr^y d [r^y - 5,7) _ d {r^/ - 6kn) 





dcki deij ' deki ddj 




'^iJi ft^, f = 1, ^1 




dr^ d mxr/A') d (r^^ - 5.,^) d {[V]xy/A) 




^0-33 deij ' 5cr33 dcij 




^,J = 1,2, 


dr^ 


d {[VFk,/F^,]xY/A') d {r^y - (5,,7) _ 8 {[VFk^/ F,^]xy /A) 


da^k 


deij ' dcr-sk deij 




hj,k = 1,2, 




diiVFks/FssUv/A') _ d{[V]xY/A') 



(74) 



(75) 



<9cr33 da^k 
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(76) 



(77) 



Eq. (1741) represents the effect of lateral deformations on the interface stress. Eq. (1751 ) describes 
the interfacial Poisson effect, in which lateral deformations of the interface produce changes in 
the "interface thickness" (excess volume per unit area). Because Eqs. (r74l)-(r76l) involve changes 
in interface area, their Lagrangian and physical forms are different. By contrast, in Eq. (1771 ) the 
derivatives are taken at a constant area and thus the Lagrangian and physical forms are identical. 
In such cases we present only the Lagrangian form of the relation. 

2. Mechanochemical relations 

Elastic deformations parallel or normal to the interface affect interface segregation. In turn, 
changes in segregation can produce changes in interface stress, interface excess volume and inter- 
face excess shears. We will present only Maxwell relations for substitutional components when the 
diffusion potentials Mki are varied. For interstitial components, the relations have a similar form 
but with the diffusion potentials replaced by the chemical potentials fii. The effect of deformations 
parallel to the interface on the interface segregation is described by the relations 

dr^^ d{[Nk]xY/A') a(rg^^-^,,-7) d{[N,]xY/A) 

dMki da, ' dMki de,j ' (78) 

I, J = 1,2; k = 2,...,K. 

Because the interface area changes, the Lagrangian and physical forms of this relation are different. 

The effect of the stress components a^i, a^2 and (133 on segregation is described by the relations 



d{[V]xYlA') d{[Nu]xY/A') 



dMki 9(T33 

d {[VFi,/F3s]xY/A') d {\N„]xyIA') 



A; = 2,...,K, (79) 



n,. n > ^ = 1.2; fc = 2,...,K (80) 

Since the derivatives are taken at a constant area, the Lagrangian and physical forms are identical. 



3. Thermomechanical relations 

Such relations describe the effects of temperature on interface stress, excess volume and excess 
shears. The relations generated by the adsorption equation would contain the excess entropy [S]xy 
which is not easily accessible. Instead, we will use the Gibbs-Helmholtz equation (fTOl) which does 
not contain [5']xy. The following Maxwell relations are obtained: 
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dT deij ' dT de^j ' (81) 



deij 


dT 

^,J = 1,2, 


dcij 


d{[V]xY/A'T) 


_dmxY/A'T') 





dT dass 

d{[VFk3/Fss]xY/A'T) _ dmxv/A'T') 



(82) 

.^ . , k = l,2. (83) 

dT dask 

where the potential \1' is given by Eq. (TtTI) . In Eq. (fSTI) . the Lagrangian and physical forms of the 
same relation are different. 

4. Thermo chemical relations 

Using the Gibbs-Helmholtz equation (iTOl) . we can evaluate the effect of temperature on in- 
terface segregation of substitutional and interstitial components. The corresponding derivatives 
involve the excess [^]xy instead of [S]xy and read 



d{[Nu]xYlA'T) dmxYl^T 



lrn2\ 



dT dMki 

d{[ni]xY/A'T) dmxY/A'T^) 



k = 2,...,K, (84) 



l = l,...,L. (85) 



dT dfii 

For substitutional components, the derivatives are taken with respect to diffusion potentials, 

whereas for interstitial components with respect to chemical potentials. The Lagrangian and phys- 
ical forms of these relations are identical. 

V. RELATION TO OTHER TYPES OF INTERFACES 
A. Incoherent solid-solid interfaces 

Incoherent solid- solid interfaces differ from coherent in two ways: 

1 . Lateral deformations of the two phases are allowed to be different and independent of each 
other, as long as they preserve the orientation of the interface plane. Accordingly, the defor- 
mation gradients of the phases must still have the upper-triangular form but need not satisfy 
Eqs. (1251) and (l26l) with equal lateral components. 
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2. Incoherent interfaces do not support static shear stresses applied parallel to the interface 
plane, responding to such stresses by sliding. 

These conditions can be reconciled with the coherent interface theory by considering only pro- 
cesses in which the shear stresses a^i and (T32 are identically zero and the lateral deformations 
of the phases remain equal. Under these constraints, the incoherency of the interface does not 
manifest itself and all equations derived for coherent interfaces are valid for incoherent ones. The 
number of independent variables reduces to {K + L + 3) and all equations are significantly simpli- 
fied. In particular, the equations no longer contain the components Fj3 of the deformation gradient 
F [Eq. (|38l) 1 and the latter need not be introduced. The shapes of the two-phase and single-phase 
regions used in the thought experiments become unimportant; only their volumes appear in the 
final equations. 

For example, the total interface free energy 'jA is given by the simplified form of Eq. (1481 ), 

K L 

-fA = [U]xY - T[S]xY - J2MNk]xY - J^/x^Nxy - [V]xyct33. (86) 

k=l 1=1 

The potentials (pm defined by Eq. (fTSl) reduce to 

K L 

(t>m = U/N - TS/N - J^MkmCk - Y,^'lCl - as^n. (87) 

k=l 1=1 

The adsorption equation becomes 

K L 

d (7A) = -[S]xYdT - ^[iV,]^yrf0fc - Y^MxYd^ll - [\/]xyc^a33 + Y. ^^^^^^^J- (^8) 

k=l 1=1 i,j=l,'i 

where we use the current state of one of the phases as the reference state of lateral strain. The 
interface stress tensor simplifies to 

where CTjj is the Cauchy stress tensor averaged over a region of volume V . 

It should be emphasized, however, that the above equations describe only some of the possible 
state variations of an incoherent two-phase systems. They do not include variations in which 
the phases undergo different lateral deformations and thus slip against each other. Due to such 
variations, an incoherent two-phase system possesses more degrees of freedom than a coherent one 
with the same number of substitutional and interstitial components. Thus, incoherent interfaces 
cannot be considered a particular case of coherent interfaces. They require a separate treatment, 
which will be presented elsewhere. 
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B. Grain boundaries 

Grain boundary is an interface between regions of the same crystalline phase with different 
lattice orientations. As other solid-solid interfaces, GBs can be coherent or incoherent. Coherent 
GBs can support not only stresses normal to the GB plane but also shear stresses parallel to it. 
When temperature^^ and/or chemical composition^^ change, some coherent GBs can change their 
structure to one that permits GB sliding. The GB becomes incoherent. 

Two different cases must be distinguished: when the grains are thermodynamically identi- 
cal and when they are not. By definition, the grains are considered thermodynamically identical 
when the phase-change equilibrium condition (|30l ) is satisfied as a mathematical identity once the 
equilibrium conditions (i)-(iv) formulated in Section ITII Bl are satisfied. In other words, the phase- 
change equilibrium condition need not be imposed as a separate equation of constraint. This can 
be the case when the grains are stress-free and uninfluenced by electric, magnetic or other fields. 
Thermodynamically identical grains can be treated as parts of the same single-phase system. On 
the other hand, in the presence of mechanical stresses or applied fields, the equilibrium thermody- 
namic states of the grains can be different. For example, when the solid is elastically anisotropic 
and the grains are subject to mechanical stresses, they either never reach equilibrium or can reach 
an equilibrium state in which their elastic strain energy densities and chemical compositions are 
different. In the latter case, the phase-change equilibrium condition (l30l) is not satisfied automat- 
ically and must be imposed as a separate constraint. Such cases should be formally treated as 
if the grains were two different phases. Accordingly, all thermodynamic equations developed in 
Sections Hn] and |IV] for phase boundaries directly apply to this case. 

There are situations when, due to crystal symmetry, the grains remain thermodynamically iden- 
tical even in the presence of certain mechanical stresses. As an example, consider a coherent sym- 
metrical tilt GB. In the unstressed state, the grains are identical and form a single-phase system. 
Due to the mirror symmetry across the boundary plane, the lateral deformations dcy (i,j = 1, 2) 
and the normal stress (733 leave the grains identical. Moreover, due to the twofold symmetry 
around the axis X2 normal to the tilt axis, the shear stress a^i also leaves the grains identical. 
Thus, when the system is subject to these deformations, it continues to be a single-phase system.— 
That Eq. (l30l) is satisfied in this case as an identity can be seen from the fact that the differences 
(f/^ — f/") , [S'^ — 5*") , f N^ — N^ j and [n^ — nf) related to grain regions containing the same 
total numbers of substitutional atoms are zero by the symmetry. The remaining terms in Eq. d30] ) 
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represent the mechanical work Wm and are given by Eq. (I3T| ). For a symmetrical tilt boundary 
Fj^ = F{^ and F33 = F^^; only the shear components F^ and F^ are different. But the term 
( F23 — F23 ) cr32 vanishes due to (132 = 0, resulting in Wm = 0. 

Generalizing this example, it can be stated that the grains remain thermodynamically identical 
during a variation of state of the system when: 

(i) The differences {W^ - t/"), (5^ - S"), (n[ - nA, (nf - nf) and {V^ - V^) remain 
zero for all grain regions with N^ = N"; 

(ii) The work Wm of the transformation of one grain to the other, given by Eq. (13T1) . remains 
identically zero. 

The term in W^ with i = 3 equals (V^ — V"') (T33 and thus vanishes. As a result, the condition 
Wm = reduces to the identity 



J2FuF22V'[Ff,-F,^)a,.^0. 

i=l,2 



(90) 



As discussed in Section IIII CI the left-hand side of this expression is the work of the shear stress 
along the transformation vector t projected on the interface plane. This identity is satisfied term 
by term when the components Fj3 are equal (F[^ = F^g) for the directions i in which the stress 
component a^i is nonzero. In fact, Eq. (|90l ) reduces to this case after an appropriate rotation of the 
coordinate axes. Note that coupled GBs^i can be equilibrated under stress as long as the relevant 
component of a^i is zero. 

Under the above conditions, the phase-change equilibrium equation (l30l) is satisfied as an iden- 
tity. Furthermore, it can be shown that when conditions (i) and (ii) are satisfied, Eqs. (|45T ) and 
(|46|) become identical to each other and only one of them should be solved simultaneously with 
Eq. (|4TI) . As a result, jA is obtained by solving a system of only two equations, giving 



K L 

-fA = [U]x-T[S]x-J2Mki[Nk]x-MN]x-J2MMx- Yl [VF^3/F,,]xas^, 

k=2 1=1 i=l,2,3 



where 



(91) 



[Z]x ■■- 



Z X 

Z" X" 



Z - Z"X/X". 



(92) 
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Index a refers to one of the grains. By specifying X, one term in Eq. (I9TI ) is eliminated. The 
coefficients \yFi^/F^^x are the excess shears (i = 1 or 2) and excess volume [y]x {i = 3) of the 
GB. 

Similarly, if conditions (i) and (ii) are satisfied, then the Gibbs-Duhem equations (l33l) and (|34]) 
for the grains become identical to each other and the adsorption equation is obtained by solving a 
system of only two equations: 

K L 

d{^A) = -[S]xdT-Y}Nk\xdMkl-[N]xd(|)l-Y^nl]xd^ll 

_ k=2 1=1 (93) 

- Y^ [VF,^/F^^]xda^i + Y^ T^jAdcj,. 

j=l,2,3 «J=1,2 

Again, one variable in Eq. 03\ is eliminated by specifying the extensive property X, which re- 
duces the number of independent differentials to (K + L + 6) . The actual number of independent 
variations is less due to the symmetry-related constraints imposed for preservation of the identity 
of the grains. In the absence of shear stresses, the last but one term in Eq. (|93] ) reduces to [V]x'^o"33. 
In this particular case, Eq. (|93T ) can be applied to both coherent and incoherent symmetrical tilt 
boundaries. We emphasize again that for Eqs. (1911 and (|93T ) to be valid, the condition F[^ = F^ 
must be satisfied for the directions i in which asj 7^ 0. 

The last term in Eq. (|93T ) contains the GB stress t^ . Assuming that grain a is the reference 
state of strain, it is straightforward to derive 



-s = J[y'Q^^]x 



\ fc=l,2 / 



A 
X 



(94) 



(a,^.V^"-(5,,a33V^"),2,j = l,2. 



AX^ 
Here, V is the bicrystal of volume, B is the displacement vector of the upper boundary of the 

bicrystal during the GB formation [cf. Fig. |3td)], aij is the volume- averaged stress tensor in the 

bicrystal, and all quantities with superscript a refer an arbitrarily chosen homogeneous region of 

grain a. In the particular case when X = A^ we have 






\ k=l,2 J 



^ [0%^- - <5,,a33V^") , ^, J = 1, 2. (95) 



AN^ 
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For this choice of X, expression (|9T| ) for 7A takes the form 

K L 

fc=2 1=1 

- as^[VU-AJ2B^(T■,.■ (96) 

i=l,2 

Eqs. (|951) and (l96l) will be used in Part 11 of this work.— 



VI. DISCUSSION AND CONCLUSIONS 

We developed a thermodynamic theory of coherent solid-solid interfaces in multicomponent 
systems under a general non-hydrostatic state of stress. All equations were derived directly from 
the First and Second Laws of thermodynamics. No Hooke's law or any other constitutive laws 
of elastic deformation were invoked. No assumptions were made regarding the interface structure 
other than the conservation of sites and elastic response to applied shear stresses. 

To circumvent the problem of undefined chemical potentials of substitutional components, we 
treat such components separately from interstitial components using diffusion potentials intro- 
duced by Larche and Cahn.""^^ Diffusion potentials in non-hydrostatic solids are well-defined 
quantities and, similar to chemical potentials, are uniform throughout an equilibrium system. Be- 
cause a system containing K substitutional components has only {K — 1) diffusion potentials [see 
Eq. do])], the requirement of their equality in coexisting phases (M^^ = M^^ k = 2, ...,K) must 
be augmented by one more condition, namely, the phase-change equilibrium equation discussed 
in Sections imi] and [mCl 

As an alternative to diffusion potentials, one can formulate the equilibrium conditions in terms 
of the 0-potentials introduced in this work [Eq. ([181)1. Equilibrium with respect to substitu- 
tional components is then expressed by K relations 0^ = (j)^ (k = 1, ...,K), which subsume 
the phase-change equilibrium condition. Written in terms of the 0-potentials, many equations of 
phase equilibrium and interface thermodynamics look similar to the familiar equations for fluid 
systems,- with the 0-potentials playing the role of chemical potentials. In the particular case of 
hydrostatically stressed solids, the 0-potentials coincide with real chemical potentials, which are 
well-defined thermodynamic quantities under hydrostatic conditions. It should be emphasized, 
however, that the 0-potentials do not solve the problem of undefined chemical potentials in non- 
hydrostatic solids in general. If we choose a different interface plane between the same two phases, 
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the 0-potentials will need to be redefined and will generally take on different numerical values. 

An important result of this paper is the coherent phase coexistence equation derived in Section 
IIII D[ It can be written as 

K L 

= -{S}xdT-Y,{Nk}xd(t>k-Y.^nl}xd^il 

fe=i 1=1 (97) 

4=1,2,3 i,i=l,2 

and is a generalization of the Clapeyron-Clausius equation to non-hydrostatically stressed mul- 
ticomponent systems. In the particular case when the phases are hydrostatic, we have Qij = 0, 
0fc = /^fc> and this equation reduces to 

K+L 

= -{SjxdT - Y {Nk}xdfik + {V}xdp, (98) 

k=l 

where we use the same symbol Nk for the amounts of substitutional and interstitial components. 
For single-component phases this equation recovers the classical Clapeyron-Clausius equation 
= —{S}NidT + {V}Nidp.^^ Eq. (|97] ) predicts a rich variety of relations between temperature, 
stress and chemical compositions of coexisting phases, which call for testing by experiments and 
simulations in the future. 

Hydrostatic phase coexistence conditions and the Clapeyron-Clausius equation have been ex- 
tensively tested by experiment and simulations. The conditions of coherent equilibrium between 
non-hydrostatic phases^— (Section fill Bl) and the phase coexistence equation (|971) derived in this 
work call for a similar verification. An analogue of Eq. (|97t for solid-fluid systems has been re- 
cently tested by atomistic simulations which implemented strongly non-hydrostatic conditions in 
the solid.'^'^ A similar analysis and simulations could be performed for solid-solid interfaces. In 
particular, it should be possible to test Eq. (1371 ) by studying the relation between variations in the 
diffusion potential and the shear stress. 

We defined the interface free energy 'jA as the reversible non-mechanical work performed 
in a thought experiment in which the interface was formed by transforming a part of a single- 
phase system into a different phase. We have shown that 'jA can be expressed as excesses of 
different thermodynamic potentials, depending on the choice of the reference properties X and Y 
in Eq. (|48] ). Two examples of such potentials are given by Eqs. (|52l ) and (|53] ). 

Despite the freedom of choice in expressing 'jA as an excess, it still remains the work of inter- 
face formation and therefore must be unique. The following comment is due in this connection. 
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The last sum in Eqs. (1471) and (1481) is the mechanical work performed by the applied stress during 
the phase transformation. This work depends on the transformation strain, which we assumed to 
be known. However, crystal symmetry may lead to multiplicity of possible transformation paths 
between the same two phases. In other words, given the same reference state there can be differ- 
ent deformation gradients F" and F'^ producing exactly the same physical states of the phases. 
Some of such alternate transformation paths can actually be implemented in experiments or atom- 
istic simulations. The situation is similar to the existence of symmetry-dictated multiple modes of 
coupled GB motion. Depending on the temperature and other factors, a moving GB can produce 
different shear deformations of the receding grain, each corresponding to a different mode of cou- 
pling. The existence of multiple coupling modes was confirmed by both simulations^^^— and 
experiments^— on symmetrical tilt grain boundaries. 

The multiplicity of possible transformation strains seems to create the following paradox: given 
a two-phase system in its current state, how does one know which of the transformation strains 
should be plugged in Eq. (|48l) to obtain the correct value of 7^4? The answer lies in the fact that the 
phase equilibrium equation (l35l) also depends on the transformation strain through the coefficients 
in front of da^i. As discussed in Section ITlI CI one of the equilibrium conditions is the neutrality 
with respect to spontaneous interface displacements in an open system. The actual deformations 
of the phases occurring during such displacements determine the equilibrium states of the phases. 
Different deformation gradients will lead to different equilibrium states. The equilibrium states of 
the phases, in tum, affect the interface free energy and all other interface properties. If the actual 
transformation strain realized in a given experiment or simulation changes, so will the equilibrium 
states of the phases and thus 'jA. Thus, knowing the current state of the two-phase system, one 
should be able to tell which of the symmetrically possible strains is actual and thus correctly 
predict the interface properties. 

Larche and Cahn^^ — realized that thermodynamic equilibrium between solid phases depends 
on the degree of coherency of the interface. The above discussion suggests that, for a coherent 
interface, the equilibrium also depends on the particular coherent transformation strain (out of 
several possible by crystal symmetry) which is realized during the interface displacements. So do 
all thermodynamic properties of the interface. Mathematically, the hypersurface of phase coexis- 
tence in the parameter space can have multiple sheets corresponding to different transformation 
strains between the same two phases. Switches between the transformation strains can cause 
abrupt changes in states of the phases and thus in interface properties, and are similar to interface 
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phase transformations. The strain multiplicity and its consequences for interface thermodynamics 
are worth exploration in future work. One way to do so would be to use Eq. (|37l) and apply the 
stress in different directions which activate different transformation paths. 

The generalized adsorption equation derived in this paper expresses the differential of the inter- 
face free energy in terms of a set of independent intensive parameters characterizing the equilib- 
rium state of a coherent two-phase system. Different forms of this equation are given by Eqs. (|58] ). 
(|59l ). (|60l) . (|65l) . (I66l) and (TTOl) . The adsorption equation can be considered as the differential form 
of the fundamental equation of the interface, whereas Eq. (|48l) the fundamental equation is a func- 
tional form. The differential coefficients in the adsorption equation define those interface excesses 
which are measurable physical quantities. In particular, the interface stress tensor r emerges from 
the coefficients in front of the lateral strains eij and is formally defined by Eq. (|6T| ). It should 
be emphasized that the excess formulation of the interface stress presented in this work is not as 
trivial as for hydrostatic systems. In the latter case, the lateral stresses in the phases are equal and 
their interface excess is unique. In non-hydrostatic systems, the lateral stresses in the phases are 
different and their excess, if calculated relative to a dividing surface, depends on its placement. In 
terms of Cahn's generalized excesses,- r depends on the choice of the reference properties X and 
Y. In this sense, the interface stress between two non-hydrostatically stresses solids in not unique. 

To understand the origin of this non-uniqueness, it is instructive to consider its Lagrangian for- 
mulation (l63l) . As indicated by Eq. (I64l) . tl is the response of the interface free energy 7^ (per unit 
reference area) to elastic deformations of the interface. Such deformations must be implemented 
in such a way that to preserve the phase equilibrium. In other words, the derivative in Eq. (|64]) 
must be taken along a certain direction on the phase coexistence hypersurface in the parameter 
space. Derivates taken along different directions give generally different values of the interface 
stress, resulting in its multiplicity. The direction in which we take the derivative is controlled by 
the choice of the reference properties X and Y. 

Another excess quantity appearing in the generalized adsorption equation is the interface excess 
shear. It characterizes the local elastic shear deformation of the interface region in response to a 
shear stress applied parallel to the interface. The excess shears [VFi3/Fs3]xy and [VF23/Fs3]xy 
are conjugate to the shear stress components a-^i and (732, respectively. Clearly, the excess shears 
are specific to coherent interfaces and are undefined for interfaces that do not support shear 
stresses. This explains why they did not appear in previous versions of the absorption equation ex- 
isting in the literature. For practical purposes, the excess shears can be normalized by the physical 
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or Lagrangian interface area, giving the specific shears [l^Fjs/FssjxyM ^^^ [VFi^/ F^^Ixy /A' , 
respectively. These quantities are similar to the GB "slip" introduced in the context of the effective 
elastic response of GBs in polycrystalline materials.— Excess shears of individual GBs in copper 
calculated by atomistic methods will be reported in Part II of this work.— 

The analysis presented in this work is limited to a plane interface between semi-infinite homo- 
geneous phases. In the future, some of these constraints could be lifted by including, for example, 
the effect of curvature or inhomogeneity. Such generalizations appear to be extremely challenging 
but could benefit from the ideas and methods developed previously in the mechanical theories of 
interfaces.— — As already mentioned, the mechanical theories consider only mechanical equilib- 
rium between the phases and do not impose the conditions of chemical or phase equilibrium. They 
consider deformations of an already existing interface and disregard the process of its formation, 
which, as we saw above and will see again in Part II,— can be accompanied by finite transfor- 
mation strains. This prevents the mechanical approach from properly defining 7 and deriving the 
adsorption equation. Nevertheless, various interface excess stresses and strains were identified, 
and carefully described mathematically, under much more general conditions than in the present 
work. Besides the excess shear and excess volume considered here, the mechanical analysis re- 
veals a few other excess strains arising, for example, when the phases are curved and/or capable 
of relative rotation and/or tilt.— 

Finally, we note that our analysis neglects the existence of vacancies, which is justified by their 
small concentration in most solids. In the absence of vacancy sources and sinks in the system, the 
total number of vacancies is conserved and they could be included in our analysis as simply one 
of the substitutional components. In this case, the number of vacancies could be treated as one 
of the independent parameters alongside the amounts of real substitutional components. As an 
alternate model, the vacancies can be assumed to be in equilibrium with some sources and sinks 
existing far away from the interface. In this case, the number of vacancies in any reference region 
of the system is a dependent parameter, whose value can be determined from the condition of 
equilibrium with the sources and sinks. It should be noted, however, that the vacancy equilibrium 
depends on specific properties of the sources and sinks. For example, one can assume that the 
phases terminate at surfaces parallel to the interface. Suppose the surfaces are in contact with an 
inert atmosphere exerting a pressure p and are capable of absorbing and creating vacancies. The 



number of vacancies in such a system is readily predictable and depends on p (Ref. 



34|) . However, 



the shear stresses 0-31 and 0-32 have to be zero because of the presence of surfaces. Thus, this model 
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will not capture the interesting interface properties associated with the shear stresses. 

On the other hand, a uniform distribution of vacancy sources and sinks inside the phases, e.g. 
in the form of climbing dislocations, would require a radical revision of the underlying assump- 
tions of the present analysis, particularly regarding the conservation of sites. This would also raise 
the questions of possible creep deformation of the stressed phases and the legitimacy of using the 
reference state formalism for the description of elastic deformations. In view of these complica- 
tions, analysis of the possibility of incorporation of equilibrium vacancies in thermodynamics of 
coherent interfaces is left for future work. 
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is not discussed here but it is straightforward to generalize our analysis to such solids. 
^^ In their analysis of coherent equilibrium, Larche and Cah n'^-^^ used the first Piola-Kirchoff tensor which 

is the transpose of the tensor used in our work. 
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stress causes coupled GB motion,— which cannot be prevented unless we create different chemical 
compositions in the grains and thus a thermodynamic driving force balancing the driving force of coupled 
motion. But then the grains essentially become two different phases. 
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Figure 1 : Two-dimensional schematic of a volume element undergoing a finite deformation. In the reference 
state (dashed lines), the volume element is a unit square. The components of the deformation gradient F 
correspond to the new lengths and projections of the edges of the reference square in the deformed state 
(sohd hnes). 
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Figure 2: Two-dimensional schematic of phases a and /3 separated by a coherent interface. When the 
interface moves down, the striped region of phase a shown in (a) transforms the striped region of phase 
/3 shown in (b). In (c), the reference volume element V' (shaded unit square) transforms to the coherent 
phases a (dashed lines) and /3 (solid lines). The differences between the deformation-gradient components 
Fi3 and F33 in the phases form the transformation vector t. 
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Figure 3: Two-dimensional schematic of coherent transformation of a region of phase a to a two-phase 
region containing phases a and /3 separated by an interface, (a) Reference state of the region, (b) Deformed 
phase a. (c) Two-phase region, (d) Overlapping shapes of phase a and the two-phase region, showing the 
displacement vector B. The open circles mark selected physical points labeled a through d with the prime 
indicating the reference state and the asterisk indicating the two-phase state. The parallelogram defined by 
the vertices a, b, c* and d* is shown separately in (e). 
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Figure 4: The two-phase region used in the derivation of expressions for the total interface free energy 7 A. 
The single-phase regions used in the derivation are chosen outside the two-phase region. 
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